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Abstract 



In relativistic potential models of quarkonia based on a Dirac-type of equation with a 
local potential there is a sharp distinction between a linear potential V which is vector-like 
and one which is scalar-like: There are normalizable solutions for a scalar-like V but not for 
a vector-like V. It is pointed out that if instead one uses an equation of the no-pair type, 
■ which is more natural from the viewpoint of field theory, this somewhat bizarre difference 

Ojl . disappears. 

in 
a: 

1. Since the discovery of the narrow resonances in the GeV region, interest in a potential model 
description of these mesons and less charming ones as quark-antiquark bound states has continued 
unabated. A long-standing problem which arises in this connection is the following: If one tries 
to include relativistic effects with a Dirac-type of equation involving a purely local potential there 
is a dramatic difference between a linear potential which is vector-like, V = k v r, and one which is 
scalar-like, V = fiifakgr. For the vector-like case there are no normalizable solutions. In view of 
the continued interest in such models|], Q], it may be useful to point out that if one uses equations 
of the no-pair type || , which are much more natural in the context of field theory, this dichotomy 
is not present. 

By way of review and for simplicity let us consider the case of an antiquark much heavier than 
the quark and take as a starting point a Dirac equation of the form 

(h D + V)i) = Eip, h D = a ■ p op + pm, (1) 

with V linear in r. A simple way to see the trouble which arises for V = k v r is to decompose the 
wave function ip into a sum 

^ = ^,+ + ^-, ^ = (3^ (2) 
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where the /5 ± are zero-momentum projection operators, defined by 

P± = (1 ± (3)/2. (3) 

In the standard representation of the Dirac matrix f3 this is essentially a decomposition into upper 
and lower components, but no use need be made of this fact. From (1) and (2) we have 

if>- = (E + m- fc„r) _1 a -pop (4) 

With E > 0, this implies that ip~ has a pole singularity at r = (E + m)/k v and is therefore not 
integrable. Since 

< >=< ifj + \ip + > + < ifj~\ip~ >, (5) 

the norm of ip~ will be infinite even if that of ip + is finite. However, if the potential is scalar-like, 
V = j3k s r, the minus sign in the denominator in (4) changes to a plus sign, 

ij~ = (E + m + k s r)- 1 a-p op ij + , (6) 

and there is no singularity. Thus, if both scalar and vector confining potentials are used it is 
necessary to have k s > k v . The same feature holds in the two-body equations of a similar type. 

2. The corresponding no-pair equation does not suffer from this dichotomy. The counterpart 
of (1) is now 

(h D + A + p UA° + p )4> + = Eij + (7) 
where A+ is the positive-energy Casimir projection operator, defined by 

A op = (E op + h D )/2E op , £ op = (p2 p + m 2)V2. (8) 
The subscript " +" indicates that ip + satisfies 

m + = ^ (9) 

and is thus a superposition of only positive-energy plane waves. From (9) it follows that 

i>+ = Rop^t, R op = {E op + m)~ l a ■ p op (10) 

Thus, regardless of the choice of U, there is now no singularity involved in the equation relating 
ip\ an d Vv Furthermore, we have 

<^\^> = <^X\Rl p RoM> 
and in p-space the operator Rl p R op is just p 2 /[E(p) + m] 2 , which is bounded by unity. Thus 

<v>;l^;> < <^X> (ii) 

and if the norm of ip^. is finite, so is that of tp+ and hence that of ip + . 
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3. To complete the argument let us compare the Schroedinger-Pauli form of the eigenvalue 
problem for the two cases of interest. As in earlier work it is convenient to introduce a new wave 
function which in p-space differs from ip~f by a slowly varying factor |3|. [|] . 

i)% = A op (j), Aop = (E op + m)/2E op . (12) 

Then ip + = ip+ + ip~ = (1 + R op )ip^, or 

V+ = S0, S = A op (l + R op )f3 + . (13) 

It is easy to verify that because of the extra factor S is pseudo-unitary, S^S = [3 + . Since 

/3+0 = (14) 

it follows that ip + and have the same norm. On multiplying (9) on the left by S> one finds that 
satisfies the equation 

H red <P = £0 (15) 

where 

H red = n S^HS n . (16) 

The quotes indicate that f3 is to be replaced by unity when acting directly on 0. For a potential 
U of the generic form 

U = U V + (3U S (17) 

computation yields 

H re d = E op + V re d, (18) 

where 

Vred = A op U+A op + {2E op A op )- 1 a ■ p op U-a ■ p op (2E op A op )- 1 (19) 

with 

U± = U v ± U a . (20) 

Since 

a ■ p op U^a ■ p op = p op • U-p op + a ■ (gradU_) x p op , 

the main difference between a pure vector and pure scalar potential is a change in sign of part 
of the spin-independent relativistic correction and in the sign of the spin-orbit interaction. Since 
these corrections do not dominate the effective interaction, one expects that there are normalizable 
solutions both in the scalar case and in the vector case [El. 



4. Of the making of potentials, as for books, there is no end. One criterion in a semi- 
phenomenological analysis of systems for which it is makes sense to attempt a description in 
terms of relativistic Schroedinger-like equations is simplicity. The use of purely local potentials 
lends itself to this because it limits the proliferation of parameters. Another criterion is to take 
note of the implications of field theory. For two spin-1/2 particles three-dimensional equations 
tied to field theory inevitably lead to effective interactions which involve projection operators. A 
reasonable compromise is therefore to consider equations of the no-pair type 0, ||: 

[(«i • Pop + Pimi) + (-a 2 • Pop + /32m 2 ) + V ++ \if) = E4> (21) 
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where 

V ++ = A ++ UA ++ , (22) 
A ++ is the projection operator product A^(l) A^(2) and 

A?(i)V = V (< = 1,2). (23) 

One may choose U to be purely local without running into difficulties. Note that the nonlocality 
of the projection operators does not introduce any new parameters, since it involves only the 
constituent masses, already present in the Dirac Hamiltonians. 

As confirmation of the fact that no problems arise even if U is purely scalar-like, it should be 
noted that no difficulties are encountered with the numerical solution of (21) when U is chosen to 
have the scalar form || 

U = kfafcr. (24) 
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